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Abstract 

We present measurements oi D —^ Kgn and D -^ Kl'^ branching fractions using 281 pb~^ of 
V'(3770) data at the CLEO-c experiment. We find that B{D° -^ Ks'^^) is larger than B{D° -^ 
KItt^), with an asymmetry of R{D^) = 0.108 ± 0.025 ± 0.024. For B{D+ -^ Ks^'^) and B{D+ -^ 
K^TT^), we observe no measurable difference; the asymmetry is R{D^) = 0.022ib0.016ib0.018. The 
D^ asymmetry is consistent with the value based on the U-spin prediction A{D^ -^ K^7r^)/A{D^ -^ 
K^TT^) = — tan^ 9c, where 9c is the Cabibbo angle. 



As the dominant decay of the charm quark is to the strange quark, the final states of 
D^ and D~^ meson decays typically include K~^, K~ , Kg, and K^ mesons. While there 
have been many measurements of D decays to final states containing K"^ and Kg, until now 
there have been no measurements of decays to final states containing a K^. Typically it 
has been assumed that the branching fraction for a decay D -^ K^X will equal that for 
the corresponding decay D -^ K^X. However, as pointed out by Bigi and Yamamoto [l|, 
interference between Cabibbo-favored transitions (producing an s quark, and thus a K^) and 
doubly- Cabibbo-suppressed transitions (producing an s quark, and thus a K^) can lead to a 
difference in the rates for D -^ K^X and D -^ KgX. Here we present first measurements of 
two D decays to i^£, D^ -^ K^tt^ and D^ ^ K^n^, and we compare the branching fractions 
with those for D^ —>■ K^tt^ and D~^ —>■ K^tt^. (Throughout, charge-conjugate modes are 
implied, except where noted.) These comparisons provide information about amplitudes 
and strong phases in D ^ Ktt decays. 

For these measurements we use a 281 pb^^ sample of e~^e~ —>■ iIj{3770) events, produced 
by the CESR-c storage ring and recorded with the CLEO-c detector. The CLEO-c detector 
is a general purpose solenoidal detector which includes a tracking system for measuring 
momentum and specific ionization (dE/dx) of charged particles, a Ring Imaging Cherenkov 
detector (RICH) to aid in particle identification, and a Csl calorimeter for detection of 
electromagnetic showers. The CLEO-c detector is described in detail elsewhere |2|, y, |j] . 

The '?/'(3770) resonance is below the threshold for DDn, and so the events of interest, 
e"'"e~ -^ '?/'(3770) -^ DD, have D mesons with energy equal to the beam energy and a unique 
momentum. Thus, for identifying D'^ and D~^ candidates, we follow Mark HI [5(] and define 
the two variables AE and beam-constrained mass Mbc by: 



i y i 

where Ei and p, are the energies and momenta of the D decay products. For true D 
candidates, AE will be consistent with zero, and Mbc will be consistent with the D mass. 

We measure the branching fractions for the decays D —>■ K^tc by directly reconstructing 
the final-state particles, where the Kg is reconstructed from Kg ■—>■ tt+tt^. The decay 
D^ — > KgTT^ is measured by a separate CLEO-c analysis [6|, but CLEO-c has not previously 
measured D° — > Kgir^. In this paper, we cite the D^ — > Kgii^ result and present a 
measurement of D^ -^ Kgir^. 

The decays D — > K^tt have not been previously measured due to the difficulty of K^ 
reconstruction. Fortunately, the clean DD environment allows us to measure these decays 
without directly detecting the K^. Instead, we reconstruct all particles in the event except 
for the Kl - that is, a tag D and a vr - and infer the presence of a K^ from the missing 
four-momentum. Our signal is a peak in the missing mass squared distribution at the K^ 
mass squared. 

The situation for D^ decays has an added complication. When D^ and D^ are pair- 
produced through a virtual photon (J^*" = 1 ), they are in a quantum coherent state. 
Therefore the decays of D^ and D^ are subject to interference. This interference has no 
effect on the overall rate for any particular D^ or D^ decay, but it does alter how often a 
particular D^ decay occurs in combination with a particular D^ decay. Therefore, when D'^ 
decays are measured with a reconstructed tag D^^he apparent "branching fractions" of the 
D^ will vary according to the decay of the D^ [7|. This effect is especially large for CP 
eigenstate modes, such as Kgii'^ and K^^tt^. 



TABLE I: Untagged (vs. X) and tagged (vs. /) efficiency-corrected yields for C = —1 D^D^ 
events, to leading order in the mixing parameters. A^ is the number of D^D^ events, Bi is the 
branching fraction for mode i for an isolated D^ , and zj = 2 cos 6f . 



X f 



S+ 2NBs+ NBfBs^ (1 + Rwsj + rjZf + y) 
S- 2NBs- NBfBs_ (1 + Rwsj - rfZf - y) 
X 2NBf{l + Rwsj) 

The quantum correlation effects are shown in Table [H where X stands for all modes 
combined, / stands for a flavored mode, 5"+ stands for a CP-even mode, S^ stands for a 
CP-odd mode, Rwsj is the wrong-sign decay ratio B{D -^ f)/B{D ^ f), y is the D^-D^ 
mixing parameter y = Ar/2r, and rfe~''^f = {f\D^) / {f\D^). An untagged measurement 
is not altered relative to a measurement using an isolated D^ . However, measurements of 
K'gTi^ (CP-odd) and K^ti^ (CP-even), tagged by a flavored Z)° decay, are altered by factors 
of (1 + Rwsj T ^fZf =F y)- These factors depend on the tag mode, and Zf = 2 cos 5/ is 
generally not known. Since D^ —>■ K^-k^ must be reconstructed with a tag Z)°, we must 
determine the factor for each tag mode, /. We do this by comparing tagged and untagged 

Our procedure is the following: We first measure the branching fraction B{D^ -^ -^5^°) 
by reconstructing this decay without tagging a Z)°. Next, we measure the "branching frac- 
tion" for D^ -^ K^T^^ ^ with three different flavor tags. Each gives us B{D^ -^ Kg7r^){l — Cf), 
where Cf = {rjZj + y)/{l + Rwsj)- Using i3(D° -^ -^5'""°) from the untagged measure- 
ment, we obtain Cf for each flavor tag. Finally, we measure the "branching fraction" for 
D'^ -^ K^j^-K^ 1 with the same three flavor tags. Each gives us B{D^ -^ K^7r'^){l + Cf). Using 
the calculated values of Cf, we obtain B{D^ -^ ^l'^'^) from each of the three tags. These 
measurements are then averaged for the final result. 

We first measure D^ — > K^tt^ without searching for a tag D^. Candidates for D^ — > Kgii^ 
are formed by combining a K^, reconstructed by a pair of charged tracks through the decay 
Kg -^ TT+TT^, and a 7r°, reconstructed from a pair of photons detected in the Csl calorimeter. 
The invariant mass of Kg and vr" candidates is required to be consistent with the known 
mass, and tt^ candidates are then constrained to the known mass. 

Both beam-constrained mass and AE are required to be within 3 standard deviations of 
the nominal value. If there are multiple candidates in one event, we accept only the one 
whose beam-constrained mass is closest to the nominal D^ mass. Two sideband subtractions 
are used to remove background. First, a AE sideband subtraction is used to remove the 
continuum and combinatoric background (a 13% effect). Then, a Kg mass sideband sub- 
traction is used to remove the background from P)° -^ 7r+7r^7r° events in which M^ti^ti^) 
happens to be within the Kg mass window (a 4% effect). The resulting yield is 7487±101 
events. This yield is divided by the detection efficiency, 29.3%, to determine the number of 
Z}° — >■ Kgir'^ events produced. The efficiency is determined from Monte Carlo simulation, 
with a correction for 7r° detection efficiency; this correction is determined by comparing vr^ 
efficiencies measured in data and in our simulation. Finally, we use the total number of 
D^D^ events in our sample, 1.031 x 10® (from a separate CLEO-c analysis J6|). Dividing the 
efficiency-corrected yield by twice this number gives the branching fraction. 

Systematic uncertainties considered include those from: the AE cut (±0.5%), the AE 



TABLE II: Efficiencies, yields, and results for tagged D^ —f Kgir^ study. No systematic uncertain- 
ties are included in the quoted results. 



Tag mode 


K+TT- 


K+vr-vrO 


K+TT-7r-7r+ 


Efficiency 


31.74% 


31.29% 


29.97% 


Tag yield - raw 


48095 


67576 


75113 


Sideband subtracted 


47440 


63913 


71040 


Signal yield - raw 


172 


248 


276 


Sideband subtracted 


155 


203 


256 


Tag bias correction 


1.000 


1.014 


1.033 


B{Kl7r^){l-Cj){%) 


1.03±0.09 


1 1.00±0.09 1.16±0.08 



sideband subtraction (±0.8%), tracking efficiency (±0.6%), K^ detection efficiency (±1.^ 
the K^ sideband subtraction (±0.3%), and the number oi D^D^ events (±1.4%). These total 
±2.5%. The largest uncertainty is due to vr" reconstruction efficiency (±3.8%). Although 
this uncertainty is large, it cancels in the computation of quantum correlation factors and 
in the comparison of the D^ -^ Kgir^ and D^ — > K^n^ branching fractions. Therefore we 
keep it separate from the other uncertainties. 

We find a branching fraction B{D° -^ i^g7r°) = (1.240 ± 0.017 ± 0.031 ± 0.047)%, where 
the last uncertainty is from the 7r° efficiency. 

Having determined B{D^ — > Kgir^), we now measure this decay with three different tag 
modes to obtain the quantum correlation factors (1 — Cf). The three tag modes we use 
are D^ — »• K~^7r~, D^ — >■ K~^n~7r^, and Z)° — > K^ti^'k^'k^. The tag Z)° is required to be 
within 3 standard deviations of the nominal values of AE and Mbc- We select at most one 
candidate per flavor per tag mode; when multiple candidates pass our requirements, we keep 
the one with Mbc closest to the nominal D^ mass. We remove fake tag Z)° candidates by 
subtracting the AE sideband of the tag. 

In the tagged sample, we reconstruct D^ —>■ K^tt^ in the same way as in the untagged 
case. To remove fake D° -^ Kgir^ candidates, we subtract a Kg mass sideband. (No AE 
sideband subtraction is necessary since, with a tag, the Kgir'^ signal is essentially free of 
combinatoric background.) 

Although to flrst order the efficiency of reconstructing the -D° tag cancels in the branching 
fraction calculation, simulations indicate a slightly larger efficiency for D^ reconstruction 
when the signal D^ decays to K^tt^. This bias stems from the lower-than-average multiplicity 
of particles in D^ —* K^tt^ events. We obtain correction factors for these small biases from 
Monte Carlo studies. 

With the efficiencies from Monte Carlo simulations and the yields in signal and sideband 
regions, we compute the branching fractions, times quantum correlation factors, in Table [TTl 

The systematic uncertainties are similar to those in the untagged measurement. Track, 
Kg, and 7r° reconstruction uncertainties are the same, and they will cancel in the ratio of 
the tagged and untagged results. The only systematic uncertainties from the tag D^ are for 
the AE sideband subtraction and the tag bias correction factor; any other discrepancies in 
the Monte Carlo simulation would have the same effect on the tag and signal yields. 

Finally, we divide these results by B{D^ -^ f^^vr^), from the untagged measurement, to 
obtain the three quantum correlation factors (1 — Cf), where / represents the tag mode. 

We measure the D -^ K^tt branching fractions with a missing mass technique. We 



reconstruct the tag Z) in 3 Z)° modes and 6 D~ modes, and we combine it with a vr*^ or vr"*" 
to form missing mass squared: M^j^g = (pevent — Pd ~ VttY- To improve resolution, the tag 
D is constrained to have the expected three-momentum magnitude. The D — ^ K^tt signal 
is a peak in M^j^g at the K^ mass squared (~ 0.25 GeV^). 

To remove D -^ K'^ir events, as well as other backgrounds, we require that the event 
contain no extra tracks or vr^'s beyond those used in the tag D and the tt. This veto removes 
about 90% oi D ^ ^s'^ events and a few percent oi D -^ K^n events. For D" — > i^^7r° 
only, we also remove an event if it contains an extra rj —>■ '-fj. This removes much of the 
D^ -^ r]7i^ background. To determine systematic uncertainties from the appearance of fake 
extra particles in signal events, we compare how often they appear in data and in our 
simulation, using events in which both D and D were fully reconstructed. 

As in the tagged D^ —^ K^tt^ study, the tag D reconstruction efficiency is higher when 
the D decays to K^tt; therefore we apply correction factors determined from Monte Carlo 
simulations. The efficiency for observing D -^ K^tt, given that the tag was found, is also 
determined in these simulations. It is essentially the efficiency for finding the tt without any 
fake extra particles. 

For the D^ — > Ki^tt^ branching fraction measurement, the same three Z)° decay modes 
are selected with the same requirements as in the tagged D^ — > K^tt^ study (except for 
a minor difference in the order of applying cuts for the K'^tt^tt^ tag, which results in a 
slight difference in number of tags). Combining these Z)° candidates with vr'' candidates and 
rejecting events with extra tracks, vr'^'s, or r^'s, we obtain the M^j^g plot shown in Fig. [H 

A number of backgrounds slip through our extra track, 7r°, and rj vetoes and appear in 
the M^jgg plot. The modes Kgir'^ and rjir^ appear as peaks at essentially the same location 
as irOyrO; 7r°7rO peaks at M^-^^ ^ 0.0 GeV^; and i^*°7r° peaks at 0.8 GeV^. Monte Carlo 
simulations of these backgrounds are shown in Fig. [1] Other, lesser backgrounds also appear 
to the right of the i^°7r° peak. 

To determine the signal and estimate the background, we define a M^j^g signal region 
0.1 to 0.5 GeV^ as well as low and high sidebands: -0.1 to 0.1 GeV^ and 0.8 to 1.2 GeV^. 
The backgrounds are split into three groups: D^ -^ K^ir^ and Z)° — > r/vr*^; D^ —>■ tt^tt^; and 
all other backgrounds. For D^ —>■ K^tt^ and D^ —>■ t^tt^, we use Monte Carlo simulation 
to determine efficiencies for the background subtraction. For D^ —>■ tt^tt^, we scale the 
contribution to the signal region according to the yield in the low sideband. For the sum 
of all other backgrounds, we follow the same procedure with the high sideband. In total, 
about 10% of the events in the signal region are background, with half coming from K^n^, 
1/10 from each of rj7r^ and tt^tt^, and 3/10 from various other decays. 

After subtracting all the backgrounds, we obtain the yields and compute branching frac- 
tions, times quantum correlation factors, in Table lllli 

Systematic uncertainties come from the effect on signal efficiency of the veto on extra 
tracks (±0.3%), the veto on extra 7r°'s (±1.6%), the veto on r^'s (±0.5%), and the uncertainty 
in the location and width of the signal peak (±1.4%). Other uncertainties come from 
the background estimate (±1.0%), AE sideband subtraction (±0.5%), and the tag bias 
correction factor (±0.2%). These total ±2.5%. As in D^ -^ i^T^ 7r°, 7r° efficiency (±3.8%) 
is the largest systematic uncertainty; it cancels in the comparison of D^ — ^ K^tt^ and 
DO -^ Kin''. 

We have determined B{D^ -^ i('£7r°)(l ± Cf) for three different flavor tags /. Using the 
values of Cf determined from the D'^ —>■ Kgir^ measurements, we calculate B{D^ —>■ K^tt'^) 
for each tag mode. Finally, we average the results and flnd BID'' -^ -^l^°) = (0.998 ± 
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FIG. 1: Missing mass squared distribution, with all tag modes combined, for D^ -^ Xi:^, after 
removing events with extra tracks, vr'^'s, or r/'s. The points with error bars are data, and the solid 
line is a Monte Carlo simulation. The dashed, colored lines represent simulations of the peaking 
backgrounds D^ — > vr'^vr'^, Kgir^, V^^: ^^'^ K*^7r^. The difference in the peak position is due to a 
minor discrepancy in our calorimeter simulation at large photon energies; the signal region, marked 
with arrows, encompasses the peak in both distributions. 



TABLE III: Efficiencies, yields, and results for D^ 
included in the quoted results. 



ivrOyrO. 



No systematic uncertainties are 



Tag mode 


K+TT- 


K+7T-7T^ 


K+TT-TT-TT+ 


Efficiency 


55.21% 


52.72% 


49.88% 


Tag yield - raw 


48095 


68000 


75113 


Sideband subtracted 


47440 


64280 


71040 


Signal yield - raw 


367.0 


414.5 


466.5 


Background subtracted 334.8 


363.1 


418.0 


Tag bias correction 


1.000 


1.037 


1.057 


B{Kln^){l + Cf){%) 


1.28±0.08 1.03±0.06 


1.12±0.06 



0.049 ± 0.030 ± 0.038)%, where the last uncertainty is from the 7r° efficiency. 

The analysis of D^ -^ K^-k^ is similar to D^ — > K^tt^, though there are a few differences. 
Since we reconstruct a 7r+ instead of a vr", the M^j^g resolution is better. Also, we do not 
need to correct for quantum correlation. The most significant difference in procedure is that 
we perform a likelihood fit for the signal and background yields instead of counting events 
in a signal region. 

We reconstruct tag D~'s in six decay modes: D^ — ^ i^+vr"??", K^tt^tt^tt'^, Kgir^ , 
Kg7r~7i^, Kg7i'~n~n~^, and K~^K~7i~' . As before, candidates must have AE consistent with 
zero. We select one candidate per charge per mode based on the best value of AE. We 
fit the Mbc distribution for each mode to determine the number of tags, and then pass all 
candidates with Mbc near the peak to be combined with 7r+ candidates. 

The M^jgg distribution, with all tag modes added together, is shown in Fig. [2l The lines 
show a fit used to determine the signal yield. The most prominent feature is the signal peak 
at the K^ mass squared (~0.25 GeV^). A number of backgrounds are also present. First, 
fake D~ candidates produce a background which is estimated from an Mbc sideband. All of 
the other backgrounds come from other D^ decays. The largest of these are D^ — >■ Kgir^ 
(dashed, green peak under the signal), tjit^ (shoulder on the right-side tail of the signal), 
71^71^ and /U+z/^ (peak on the left of the plot), K^tt'^tt^, and vr+vr^vr^. The shapes and 
efficiencies of these backgrounds are determined from Monte Carlo simulations. The yields 
of the signal peak and the r^vr"*", tt^tt"'", and /i^t'/i backgrounds are allowed to vary in the fit; 
all other yields are fixed based on the efficiencies. 

Although Fig. [2] shows all tag modes together, we actually fit each tag mode separately. 
We calculate a branching fraction from each tag mode using the tag bias correction factor, 
efficiency, tag D~ yield, and signal D~^ -^ K^n^ yield for that mode. The tag bias correction 
varies from 1.005 (for K^t{~t{~) to 1.047 (for ii'°7r~7r~7r^). The efficiency averages to 81.6%, 
and depends little on tag mode. There are a total of 165 x 10^ tags, and a total D^ — > KIti^ 
yield of 2023±54 events. The values of the branching fraction calculated from each tag mode 
are averaged to produce the final result. 

Systematic uncertainties include those from: pion reconstruction efficiency (±0.3%) and 
particle identification (±0.25%), tag bias correction factor (±0.2%), charge of the tag D 
(±0.5%), extra track and extra vr" vetoes (±1.1%), signal peak shape (±0.7%), signal peak 
width (±1.6%), contribution of fake D^ tags (±0.4%), and M^j^^ background yields (±0.8% 
from statistical uncertainty in K^gT^^ background, ±0.3% from B{D~^ — > i^^vr"*"), and ±0.5% 
from all other backgrounds). The total systematic uncertainty is ±2.4%. 

We find a branching fraction B{D+ -^ KIti+) = (1.460 ± 0.040 ± 0.035 ± 0.005)%. The 
final uncertainty is due to the input value of B{D^ -^ K^-k^). 

To compare D -^ Kgir and D — > K^tt, we compute the asymmetries 

^ BjD ^ Kin) - BjD -. Klix) 
^ '~ B{D^ K^ir) + B{D -^ KI-k) ' 

The D° asymmetry (in which the systematic uncertainty for 7r° efficiency cancels) is R{D^) = 
0.108 ± 0.025 ± 0.024. Using B{D+ -^ i^s^+) = (1.526 ± 0.022 ± 0.038)% |Q], the D+ 
asymmetry is R{D+) = 0.022 ± 0.016 ± 0.018. 

The asymmetry between D^ -^ Kgir'^ and D'^ -^ KItt^ is consistent with SU(3) symmetry, 
and in particular the U-spin subgroup of SU(3). U-spin predicts A(Z)° -^ K^7i'^)/A{D'^ -^ 
K'^TT^) = — tan^^c", where 6c is the Cabibbo angle. This prediction is relatively insensitive 
to SU(3) breaking [81]. The amplitude ratio can also be predicted from diagrams for these 
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FIG. 2: Missing mass squared distribution, with all tag modes combined, for D^ -^ X'k^ 
removing events with extra tracks or tt^'s. The solid line shows a fit for the -D"*" 



, after 
yield. 

The many dashed lines represent the various components of the fit, added cumulatively. The small 
peak under the signal is the contribution of D^ — > K\t:^ events that are not removed by the extra 



Kl^^ 



track and vr^ vetoes. 



two processes; both have spectator and exchange diagrams which differ only by a factor of 
— tan^ Qc- However derived, the amplitude ratio implies that the asymmetry is R{D^) = 
2tan^ 9c. Using tan^c = 0.233 ± 0.001 Q, we calculate R{D^) = 0.109 ± 0.001, in good 
agreement with our measurement. 

There is no corresponding U-spin argument for the D^ decays, so no simple prediction 
is possible. Diagrams for the Cabibbo-favored and doubly-suppressed decays are differ- 
ent. Both internal and external spectator diagrams contribute to D^ -^ K^n^, while 
D^ -^ K^7c~^ has internal spectator and annihilation diagrams. Approximate predictions 
are, however, possible under certain assumptions. One analysis [lO|, based on flavor SU(3) 
with an estimate of symmetry-breaking effects, finds R{D~^) p^ 0.04, consistent with our 
measurement. This analysis also points out that the small asymmetry found for D~^ decays 
can be interpreted as a large strong phase between two contributing amplitudes in the case 
of D~^ decays, while the larger asymmetry in the D^ decays is consistent with a small strong 
phase. 
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